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Resolving quantum many-body problems represents one of the greatest challenges in physics and 
physical chemistry, due to the prohibitively large computational resources that would be required 
by using classical computers. A solution has been foreseen by directly simulating the time evolution 
through sequences of quantum gates applied to arrays of qubits, i.e. by implementing a digital quan¬ 
tum simulator. Superconducting circuits and resonators are emerging as an extremely-promising 
platform for quantum computation architectures, but a digital quantum simulator proposal that 
is straightforwardly scalable, universal, and realizable with state-of-the-art technology is presently 
lacking. Here we propose a viable scheme to implement a universal quantum simulator with hybrid 
spin-photon qubits in an array of superconducting resonators, which is intrinsically scalable and 
allows for local control. As representative examples we consider the transverse-field Ising model, a 
spin-1 Hamiltonian, and the two-dimensional Hubbard model; for these, we numerically simulate 
the scheme by including the main sources of decoherence. In addition, we show how to circumvent 
the potentially harmful effects of inhomogeneous broadening of the spin systems. 


I. INTRODUCTION 

There is a large number of problems that are well 
known to be hardly tractable with standard computa¬ 
tional approaches and resources, mainly due to the many- 
body nature of strongly correlated many particle systems. 
To overcome this limitation, the idea of a quantum sim¬ 
ulator was originally proposed by Feynman [l|: any ar¬ 
bitrary complex quantum systems could in fact be simu¬ 
lated by another quantum system mimicking its dynam¬ 
ical evolution, but under the experimenter control. This 
idea was later refined and mathematically formalized in 
quantum information perspectives by Lloyd 
Over the past twenty years, different approaches have 
been proposed to realize quantum simulators of the most 
relevant models in condensed matter physics, quantum 
field theories, and quantum chemistry Q. Most efficient 
protocols have been proposed and experimentally real¬ 
ized with trapped ions [4j, [5|. Generally speaking, quan¬ 
tum simulators can be broadly classified into two main 
categories: in digital simulators the state of the target 
system is encoded in qubits and its Trotter-decomposed 
time evolution is implemented by a sequence of elemen¬ 
tary quantum gates @], whereas in analog simulators a 
certain quantum system directly emulates another one. 
Digital architectures are usually able to simulate broad 
classes of Hamiltonians, whereas analog ones are re¬ 
stricted to specific target Hamiltonians. For a recent 
review on these different approaches, we refer to [3] and 
references therein. 

Lately, superconducting circuits and resonators have 
emerged as an extremely promising platform for quantum 
information and quantum simulation architectures 
The first and unique proposal for a general-purpose dig¬ 
ital simulator has been put forward only very recently 
0. In this proposal qubits encoded in transmons are 
dispersivcly coupled through a photon mode of a single 


resonator, and such coupling is externally tuned by con¬ 
trolling the transmon energies. However, the reported 
fidelities and the intrinsic serial nature of this setup (i.e., 
the need of addressing each pair of qubits sequentially), 
may hinder the scalability to a sizeable number of qubits. 
In addition, superconducting units are not ideal for en¬ 
coding qubits owing to their relati vely short coherence 
times. Indeed, spin-ensembles 110-1121] or even photons 
BUI have been proposed as memories to temporarily 
store the state of superconducting computational qubits. 
Here we consider an array of superconducting resonators 
as the main technological platform, on which hybrid spin- 
photon qubits are defined by introducing strongly cou¬ 
pled spin ensembles (SEs) in each resonator 0, 0. 
One- and two-qubit quantum gates can be implemented 
by individually and independently tuning the resonators 
modes through external magnetic fields. This setup can 
realize a universal digital quantum simulator, whose seal- 
ability to an arbitrary large array is naturally fulfilled by 
the inherent definition of the single qubits, represented 
by each coupled SE-resonator device. The possibility to 
perform a large number of two-qubit gates in parallel 
makes the manipulation of such large arrays much faster 
than in a serial implementation, thus making the simula¬ 
tion of complex target Hamiltonians possible in practice. 
A key novelty of the present proposal is that ensem¬ 
bles of effective S = 1 spins are used in the hybrid en¬ 
coding, which allows to exploit the mobility of photons 
across different resonators to perform two-qubit gates 
between physically distant qubits. This is done much 
more efficiently than by the straightforward approach of 
moving the states of the two qubits close to each other 
by sequences of SWAP gates, and makes the class of 
Hamiltonians which can be realistically addressed much 
larger. Long-distance operations arise whenever map¬ 
ping the target system of the simulation onto the regis¬ 
ter implies two-body terms between distant qubits. Be- 
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sides the obvious case of Hamiltonians with long-range 
interactions, this occurs with any two-dimensional model 
mapped onto a linear register, or with models containing 
./V-body terms, including the many-spin terms which im¬ 
plement the antisymmetric nature of fermion wavefunc- 
tions. 

The time evolution of a generic Hamiltonian is decom¬ 
posed into a sequence of local unitary operators, which 
can be implemented by means of elementary single- and 
two-qubits gates. Then we combine the elementary gates 
of our setup in order to mimic the dynamics of spin and 
Hubbard-like Hamiltonians for fermions. We explicitly 
report our results for the digital quantum simulation of 
the transverse-field Ising model on 3 qubits, the tunnel¬ 
ing dynamics of a spin one in a rhombic crystal field and 
the Hubbard Hamiltonian. The robustness of the scheme 
is demonstrated by including the effects of decoherence 
in a master equation formalism. Finally, we discuss the 
main sources of errors in the present simulations and the 
possibility to overcome them. In particular, we show how 
potentially harmful effects of inhomogeneous broadening 
of the spin ensemble are circumvented by operating the 
scheme in a cavity-protected regime. 


II. A SCALABLE ARCHITECTURE FOR 
QUANTUM SIMULATION 

The proposed quantum simulator is schematically 
shown in Fig. [I] It consists of a one- or two-dimensional 
(ID or 2D) lattice of superconducting resonators where 
hybrid spin-photon qubits are defined. We notice that 
large arrays of such resonators have already been shown 
experimentally In this schematic implementa¬ 

tion, qubits are encoded within square boxes. Each box 
represents a coplanar resonator containing an ensemble 
of (effective) S = 1 spins, whose collective excitations 
correspond to the transitions from the m = 0 single¬ 
spin ground state to the m = ±1 excited states, and 
can be modeled by two independent harmonic oscillators. 
Red lines represent the transition energies (continuous 
m = — 1, dashed m = 1 transitions, respectively), while 
the blue line indicates the resonator frequency in the 
idle configuration. This can be varied within a nanosec¬ 
ond time-scale by means of SQUID devices properly con¬ 
nected to the resonator [18l - l20l] . in order to match the 
spin transition frequencies. In the hybrid qubit encod¬ 
ing, a dual-rail representation of the logical units is in¬ 
troduced where the |0) M and |l) p states of qubit ^ are de¬ 
fined in the single-excitation subspace of each resonator. 
The logical state |0) p (|l) p ) corresponds to zero (one) 
photons and a single (zero) quantum in the m = — 1 os¬ 
cillator in cavity /x. This encoding has been introduced 
in previous works EMU, and it is detailed in App. A for 
completeness. The m = 1 oscillator represents an auxil¬ 
iary degree of freedom that is exploited to store the pho¬ 
tonic component of the qubit, if needed (e.g., to perform 
two-qubit gates between distant qubits, see App. B). 


(a) 



FIG. 1: (a) Elementary unit of the scalable setup, consisting 
of an auxiliary and a logical resonator. The latter includes an 
ensemble of S = 1 spins, placed at the antinodes of the mag¬ 
netic field (rotational lines) of the cavity mode. The auxiliary 
resonator contains a nonlinear element (transmon) coupled 
to the electric field of the fundamental mode, (b) Detailed 
sequence of time steps required to produce controlled-ip two- 
qubit gate between qubits /x = 2 and /x = 3 (see Appendix B 
for details). Logical cavities are represented by square boxes, 
whereas auxiliary resonators are depicted as circular boxes. 
Blue lines represent photon frequencies in the idle configu¬ 
ration (cu^(O) in the logical and <1^(0) in the auxiliary cavi¬ 
ties). The transmon (Hoi an d SI 12 ) and spin (w_i, contin¬ 
uous, and ui 1 , dashed) transition energies are indicated by 
red lines. (I) qubits are initially into state 11213 ), with the 
excitations (red arrows) stored into the photonic degrees of 
freedom (blue lines); (II) logical cavity 3 is brought into reso¬ 
nance with the auxiliary resonator j = 2, thus (III) bringing 
the photon to the auxiliary cavity. In the meantime auxiliary 
resonator 3 is detuned from the others to avoid unwanted pho¬ 
ton hoppings. In (IV) the photon is absorbed by the transmon 
(1*00,2=2) —> |gbi,y= 2 ) transition). The same hopping process 
(V) is repeated for the photon originally in cavity 2, which 
is brought to the auxiliary resonator (VI) and then absorbed 
and emitted by the transmon (| 1 / 11 , 2 ) —>• |i/> 2 , 2 ) transition) in 
a semi-resonant Rabi process (VII). The procedure is then 
repeated to bring photons back to logical cavities 2 and 3, 
leading the state back to 11213 ) with an additional phase 1 p 
acquired during the semi-resonant process. 
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The basic unit of the scalable array is represented by 
a pair of qubits connected by an interposed auxiliary 
resonator containing a superconducting transmon device 
(circular box), which is employed to perform two-qubit 
gates. It should be emphasized that this nonlinear super¬ 
conducting element is not used to encode any informa¬ 
tion, and it is left in its ground state always except during 
the implementation of the two-qubit gates. Consequently, 
its possibly short coherence times do not significantly af¬ 
fect the quantum simulation. 

In the following, we shall refer to the square boxes as 
the logical cavities labelled with Greek letters, while the 
circular ones are the auxiliary cavities labeled by Latin 
letters. Photon hopping between neighboring resonators 
is allowed by capacitive coupling. Formally, such a com¬ 
plex system can be described by the total Hamiltonian 

H = H sp i n + Htr + H p h + Hint + Hph-ph ■ ( 1 ) 

The first term describes the SEs as independent harmonic 
oscillators E (h = 1): 

Hspin — EE ( 2 ) 

m ft 

where b ^ creates a spin excitation in level m = ±1 of 
resonator p. The transmons are treated as effective three- 
level systems, with transition energies Hoi and H 12 , and 
described by 

Htr = ^2/ ^Ol|V’l,j)(V’l,il + (^12 + Hoi)|V , 2j)(V , 2j |- (3) 

j 

The time-dependent photonic term is entirely responsible 
for the manipulation of the qubits. It can be expressed 
as: 

Hph = "22 w c {t)a\dp + ^2 , (4) 

n it 

where ujg(t) = ut£(0 ) + (t) and a similar expression 

holds for <D/(t). a * (dp) creates (destroys) a single 
photon in the logical resonator p. while a| (dj) cre¬ 
ates (destroys) a single photon in the auxiliary cav¬ 
ity j. Hereafter, we will use the interaction picture, 
with Hq = H sp in + H tr + H p h(t = 0 ). Hence, within 
the rotating-wave approximation the spin-photon and 
transmon-photon coupling Hamiltonian takes the form: 

Hint = G 01 J2 [s t |V' 0 J )(V’i,i|e i( ^- noi)t + h.c. 

3 

+ + h.c. 

+ £ 

m= 1,-1 fL 

Here, the coupling constants G m for the SE are enhanced 
with respect to their single-spin counterparts by a factor 


VAf, J\T being the number of spins in the SE [22j . 
Finally, the last term in Eq. © describes the photon¬ 
hopping processes induced by the capacitive coupling be¬ 
tween the modes in neighboring cavities E2: 

H ph _ ph = -kJ2 “I+ h.c.. (6) 

(mu') 

Single- and two-qubit gates are efficiently implemented 
by tuning individual resonator modes, as shown in pre¬ 
vious works E E- Arbitrary single-qubit rotations 
within the Bloch sphere as well as controlled-phase (C tp) 
gates can be realized (see App. B for a summary). 

The present setup offers two remarkable benefits: the 
first is that using the hybrid encoding with an ensemble 
of effective 5=1 spins ensures the possibility of imple¬ 
menting Controlled-phase gates between distant qubits, 
with no need of performing highly demanding and error- 
prone sequences of SWAP gates. This is done by bringing 
the photon components of the two qubits into neighbor¬ 
ing logical resonators by a series of hopping processes (see 
App. B for details). Transferring the photons with no cor¬ 
ruption and without perturbing the qubits encoded in the 
interposed logical cavities is made possible by temporarily 
storing the photon component of these interposed qubits 
into the m = 1 spin oscillator. 

In addition, quantum simulations can be performed in 
parallel to a large degree, with resulting reduction of sim¬ 
ulation times. This is made possible by the definitions 
of the single qubits, represented by each coupled SE- 
resonator device, and by the local control of each logical 
or auxiliary resonator. Non-overlapping parts of the reg¬ 
ister can then be manipulated in parallel. For instance, in 
simulating a Heisenberg chain of N spins s = 1/2, the N 
two-qubits evolutions which appear at each time-step in 
the Trotter decomposition are performed simultaneously 
first on all N/2 ’’even” bonds and then simultaneously 
on the remaining N/2 ’’odd” bonds. Thus the simula¬ 
tion time of each Trotter step does not increase with N. 


III. NUMERICAL EXPERIMENTS 

While it is obvious that a universal quantum computer 
can be used in principle to simulate any Hamiltonian, the 
actual feasibility of such simulations needs to be quan¬ 
titatively assessed by testing whether the complex se¬ 
quences of gates needed are robust with respect to errors 
due to decoherence. Here we numerically solve the den¬ 
sity matrix master equation for the model in Eq. © with 
the inclusion of the main decoherence processes, i.e., pho¬ 
ton loss and dephasing of the transmons El ( see App. C 
for details). The potentially harmful effects of inhomoge¬ 
neous broadening in the SE will be addressed in the next 
section. 

In the following, we will consider the fidelity 
T = x/{if\p\if ), 


a'pbm^-^ + h.c. 


Hi 


.(5) 


( 7 ) 
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as a valuable figure of merit for the target Hamiltonians 
to be simulated, where p is the final density matrix and 
| ip) the target state. For the simulations shown in the 
following, we have chosen these operational parameters: 
wi/2tt = 37 GHz, w_i/2t r = 35 GHz, w c (0)/2tt = 31 
GHz, £ c (0)/2tt = 28 GHz and H 0 i/2tt = 21.7 GHz, 
Hi2/27r = 19.6 GHz (see the level scheme inside each cav¬ 
ity in Fig. [T]). We also assume realistic values of the SE- 
resonator G ±i = 40 MHz, transmon-resonator Goi = 30 
MHz, G 12 = 40 MHz and photon-photon n = 30 MHz 
couplings, respectively I'd [23]. The transmon parame¬ 
ters correspond to a ratio between Josephson and charge 
energies Ej/Eq = 25 J24|. In this regime the dephas¬ 
ing time Tf exceeds several ys while keeping a 10% an- 
harmonicity. The two chosen spin gaps can easily be 
achieved with several diluted magnetic ions possessing a 
S' > 1/2 ground multiplet, just by applying a small mag¬ 
netic field along a properly chosen direction. We have 
chosen resonator frequencies ui c and cJ c larger than usual 
experiments (e.g., twice the typical frequencies reported 
in Ref. [23h . since this helps improving the maximal fi¬ 
delity of gates. However, we emphasize that the results 
do not qualitatively depend on these specific numbers. 
For instance, in the next Section we show that remark¬ 
ably good fidelities can be obtained by using realistic and 
state-of-art experimental frequencies ( 23) or smaller. 


A. Digital simulation of spin Hamiltonians 

Since most Hamiltonians of physical interest can be 
written as the sum of L local terms, our quantum com¬ 
puting architecture can be employed to efficiently sim¬ 
ulate the time-evolution induced by any target Hamil¬ 
tonian of the type TL = 'Yln/'R-k- The system dynamics 
can be approximated by a sequence of unitary operators 
according to the Trotter-Suzuki formula (h = 1): 

U{t) = e~ iiLt m (e- ililT ■ ■ ■ e" i?iir ) n , (8) 

where r = t/n and the total digital error of the present 
approximation can be made as small as desired by choos¬ 
ing n sufficiently large jj2] . In this way the simulation re¬ 
duces to the sequential implementation of local unitaries, 
each one corresponding to a small time interval t/n. The 
set of local unitary operators can be implemented by a 
proper sequence of single- and two-qubit gates. 

The mapping of s = 1/2 models onto an array of qubits 
is straightforward. We consider here two kinds of signifi¬ 
cant local terms in the target Hamiltonian, namely one- 
(7ia^) and two-body (R^g) t erms J with a,/3 = x,y,z. 
The unitary time evolution corresponding to one-body 
terms Ra^ = bs a is directly implemented by single-qubit 
rotations R a (br). Conversely, two-body terms describe a 
generic spin-spin interaction of the form R a / = Asi a S 2 / 3 , 
for any choice of a, /3 = x, y , z. The evolution operator, 



time 

Ti 

Td 

T {CP) 

D 

0/(1) 

Tlx 

6.4 

ns 

99.99 

% 

99.94 

% 

99.76 

% 

ui 1) 

0.5 

ns 

99.99 

% 

99.98 

% 

99.86 

% 

o/( 2 ) 

n -yy 

85.8 

ns 

99.87 

% 

99.24 

% 

98.69 

% 

o/(^) 

Ttzz 

61 

ns 

99.91 

% 

99.45 

% 

99.00 

% 

7V (2) 

Ttyz 

85.8 

ns 

99.79 

% 

99.13 

% 

98.57 

% 


TABLE I: Simulation of the elementary terms of the 
Hamiltonian. Fidelity and time required to simulate the el¬ 
ementary terms of the Spin Hamiltonian. The fidelity has 
been calculated by assuming a random initial state. The 
second and third column show a comparison between the 
ideal fidelity (calculated in the absence of decoherence) and 
the real one (calculated assuming a Lindblad dynamics, with 
Q = 10 6 and T| T = 10 ys). The implemented evolution is 
U = exp , with for = Ar = 7r/2. The last col¬ 

umn reports the fidelities with a setup operating in a cavity- 
protected regime, with Q = 10 6 and T| r = 10 ys (see Section 
IV). 

3 r , can be decomposed as (2f| 

e -i\s la s 2 flT _ ^ ^ £ 2/3 ] e -ikr ^ ( (9) 

with A = Xsizhz, u x = Ry{ 7r/2), u y = R x (3n/2), 
ii z = I. The Ising evolution operator, e ~ lXsizS2zT , can 
be obtained starting from the two-qubit Cgate and ex¬ 
ploiting the identity (apart from an overall phase) 

e -iXs lz s 2z T = [$ 1 (_ ¥ ,/2)®/ 2 ]j7 Cv [/i®$ 2 (- ¥ ,/2)] ) (10) 

where ip = At. Here $(y>) is a phase gate (see App. B). 
The time required and the fidelity for the simulation of 
each term of a generic spin Hamiltonian are calculated 
by using a Lindblad master equation formalism and are 
listed in Tabic Q] We notice that the predicted fidelities 
are very high, even after the inclusion of realistic values 
for the main decohcrence channels, especially for the pho¬ 
ton loss rate which is related to the resonators quality 
factor (Q) by r p = w^/Q. These elementary steps can 
be used to simulate non-trivial multi-spin models. 

As a prototypical example we report the digital quan¬ 
tum simulation of the transverse field Ising model (TIM) 
on a chain of 3 qubits: 

RtIM = A (s\ z S2z + §2zS3z ) + b (3l x + §2x + S3x) , (H) 

where §i a are spin-1/2 operators. Figure [2] shows the 
oscillations of the magnetization, Tr[p(si z + S 2 z + S 3 Z )], 
for a spin system initialized in a ferromagnetic configu¬ 
ration. Here p is the three-qubit density matrix obtained 
at the end of the n = 10 Trotter steps of the simula¬ 
tion. The exact Trotter evolution (continuous line) is 
compared to the simulated one (points). In particular, 
red circles represent the ideal evolution, without includ¬ 
ing any source of decoherence. Errors are, in that case, 
only due to a non-ideal implementation of the quantum 
gates (see discussion below). Conversely, green and black 
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FIG. 2: Oscillations of the magnetization in the transverse- 
field Ising model. The simulation is performed on a chain of 3 
qubits, in the case b = A/2. The plot reports the expectation 
value of the total magnetization (S z ) = tr[p(si z + S 2 z + « 3 z)] 
on the final state of the system, simulated for different values 
of the resonator quality factor (points) and compared with 
the exact evolution (line). 

circles are calculated including the most important deco¬ 
herence channels, namely photon loss (timescale 1/T M ) 
and pure dephasing of the transmon (timescale Tff ). It 
turns out that photon loss is the most important envi¬ 
ronmental source of errors [l(|, while T| r « 10 fxs [2fJ 
is sufficient to obtain high fidelities at the end of the 
simulation. Indeed, the transmon is only excited during 
the implementation of two-qubit gates. The simulation 
has been performed for different values of the resonators 
quality factor. By decreasing Q the average fidelity de¬ 
creases from 96.5% (infinite Q) to 94.6% (Q = 10 7 ) and 
84.6% {Q = 10 6 ). For high but realistic [27] values of 
Q = 10' the calculated points are close to the ones ob¬ 
tained in the ideal case (with infinite Q): in that case 
the gating errors still dominate the dynamics. Finally, 
by exploiting the auxiliary m = 1 oscillator to store the 
photon component of the hybrid qubits when these are 
idle, the effects of photon loss are reduced and the fidelity 
significantly increases. The improvement is apparent in 
Fig. [2] by comparing black circular and square points; 
the final fidelity raises from 84.6% to 92% thanks to this 
storage. We stress again that the simulation time of each 
Trotter step does not increase for larger systems contain¬ 
ing more than 3 spins. Indeed, even if more gates are 
needed, these can be applied in parallel to the whole ar¬ 
ray, independently of the system size. 

The simulation of Hamiltonians involving S > 1/2 spin 
ensembles can be performed by encoding the state of each 
spin-S onto that of 2 S qubits. As an explicit example, we 
consider a chain of S = 1 spins, labelled S,;, with nearest- 
neighbor exchange interactions and single-spin crystal- 
field anisotropy, described by the Hamiltonian 

n sl = Y, AS< ■ Si+r + [DSl + E (§*, - S* )] , 



FIG. 3: Tunneling of the magnetization in a spin-1 system. 
Line: exact time evolution of (S z ) for a single S = 1 spin with 
\D/E\ = 12, after Eq. (1121) . As it is well known, the system 
oscillates between states with opposite magnetization due to 
quantum tunneling across the anisotropy barrier. Points: dig¬ 
ital quantum simulation obtained by the time evolution of two 
hybrid qubits for different values of the resonator quality fac¬ 
tor, Q, and of the transmon coherence time, T/ r , respectively. 

which reduces to the paradigmatic Haldane case for D = 
E = 0. By rewriting each spin-1 operator as the sum 
of two spin-1/2 ones (Si a = SiAa + SiBa), 7£ s i can be 

mapped onto as = 1/2 Hamiltonian, TL s i, with twice 
the number of spins. Indeed, if each A-B pair of qubits 
is initialized into a state with total spin equal to one, 

the dynamics of 7i s i coincides with that of T~L s i and can 
be simulated along the lines traced above. A proof-of- 
principle experiment, which could be implemented even 
by the non-scalable single-resonator setup described in 
Ref. [H, would be the simulation of a single spin S = 
1 experiencing tunneling of the magnetization. In this 
simple case we find (apart from a constant term): 

'Hsl = 2 Ds z aSzB + 2 E ( S x aSxB — SyA§ y B ) ■ (12) 

Figure [3] reports the comparison between the exact and 
the simulated evolution of the magnetization, assuming 
D < 0 and \D/E\ = 12, for different values of Q and 
T 2 tr . Interestingly, quantum oscillations of (S z ) are well 
captured by the simulation even for Q = 10 5 , and the 
fidelity is practically unaffected by a reduction of trans¬ 
mon coherence time to T| r = 1 fis. 

The simulation of many-spin models with S > 1 typically 
requires two-qubit gates involving non-nearest-neighbor 
qubits. These can be handled with no need of SWAP 
gates as outlined in App. B. 

B. Digital simulation of Fermi-Hubbard models 

The numerical simulation of many-body fermionic sys¬ 
tems is a notoriously difficult problem in theoretical con- 
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densed matter. In particular, quantum Monte Carlo al¬ 
gorithms usually fail due to the so-called sign-problem 
|28(. Our digital quantum simulator setup enables to ef¬ 
ficiently compute the quantum dynamics of interacting 
fermions, even on an arbitrary two-dimensional lattice. 
Although we focus on the paradigmatic Fermi-Hubbard 
Hamiltonian, the proposed scheme can be generalized to 
the quantum simulation of several other fermionic mod¬ 
els, such as the Anderson impurity model. 

The target Hamiltonian describing a two-dimensional 
TV x M lattice of Wannier orbitals is 

tiHub = — A + U (13) 

where (fi, v) are nearest neighbors (v = /x± 1, v = fi±M) 
and c^ a are fermionic operators. In order to simulate 
this Hamiltonian with our setup, we exploit the Jordan- 
Wigner transformation to map fermion operators c^ onto 
spin ones (29l - l3l| . However if such a transformation 
is applied to the Hubbard model (1131) in more than one 
dimension, the hopping (first) term results into XY spin 
couplings whose sign depends on the parity of the number 
of occupied states that are between u and v in the chosen 
ordering of the Wannier orbitals |32[ • This aspect makes 
the simulation of a fermionic system much more demand¬ 
ing than any typical spin system, because the resulting 
effective spin Hamiltonian contains many-spin terms. To 
illustrate how we address this key issue, here we consider 
the simpler case of the hopping of spinless fermions on 
a lattice (the general case of interacting spin fermions 
is discussed in App. D). The target Hamiltonian can be 
mapped into the following spin model: 

^ = -A^(-l)^+h.c, (14) 

where a = Y ^=\+1 = E^+i^t + F We simu " 

late this n-body interaction by takin g ca re of the state- 
dependent phase, similarly to Refs. HH, [H. The sign 
factor in m is obtained by performing a conditional 
evolution of the qubits interposed between the specifi¬ 
cally addressed sites, p and 17 depending on the state 
of / 1 . This corresponds to a series of controlled-Z (CZ) 
gates between qubit p and each of the qubits 7 inter¬ 
posed between fi and v. Hence, the sequence of gates to 
be implemented at each Trotter step is the following: 

n n tkz„ i7 . ( 15 ) 

7O /^<7<i/ 

For instance, in Fig. @] we show the quantum cir- 
cuit for the implementation of i-L^f = —A(c|c 5 + c' 5 di): 
controlled-phase gates (with ip = n) between qubit \i/ji) 
and each of the qubits interposed between \ 1 jj 1 ) and | ^ 5 ), 
namely l^), l^) and 1 ^ 4 ), are sequentially performed 
before and after the central block (dashed boxes), which 
implements the XY evolution: Uxy = exp{— i<p{s\ x s^ x + 
siyS 5 y )}. The latter consists of two controlled-^ gates 


(with Lp = 2At), preceded and followed by proper single¬ 
qubit rotations, implementing respectively s x s x and s y s y 
terms of the interaction, as schematically explained in 
Fig. [I] By exploiting the high mobility of the photons 
entering into the hybrid encoding, Hamiltonian terms 
involving distant qubits can be simulated straightfor¬ 
wardly. In fact, this is a remarkable advantage with 
respect to alternative solid-state platforms for quantum 
information processing. We stress that, in spite of the 
increment in the number of gates required to address the 
sign issue, a large number of hopping terms can still be 
implemented in parallel. 


IV. QUANTUM SIMULATIONS WITH 
INHOMOGENOUSLY BROADENED SPIN 
ENSEMBLES 

In this Section we show how the present scheme can 
be made robust against inhomogenous broadening (IB) 
of the SE, which is probably the major shortcoming in 
its use to encode quantum information. A certain degree 
spin inhomogeneity is unavoidable in real SEs, and may 
result from slightly disordered spin environments or by 
random magnetic fields produced by surrounding nuclear 
magnetic moments. Due to IB, the collective (’’super¬ 
radiant” ) spin-excitation that couples to the photon held 
spontaneously decays into a quasi-continuum of decou¬ 
pled, “dark” spin modes within a timescale of order h/c r, 
a being the width of the distribution of gaps in the SE. 
A possible way to deal with IB is to revert the associated 
Hamiltonian evolution by echo techniques, e.g. by using 
external magnetic held pulses resonant with the spin gaps 
HU [ 35 ]]. While possible in principle, implementing such 
a solution within our simulation scheme would be very 
demanding. Pulses should act with very high hdelity, 
and should be controlled independently for each logical 
resonator with the proper timing to restore qubits before 
they undergo gates. Here we prefer to adopt a different 
strategy, which does not require additional resources and 
works well with the present scheme. The idea is to run 
the simulation by keeping the SE in a “cavity protection” 
regime fEI - IEl . Indeed, a strong spin-resonator coupling 
provides a protection mechanism, by inducing an energy 
gap between the computational (super-radiant) and the 
non-computational (dark) modes [38[, thus effectively de¬ 
coupling these two. This mechanism has been experimen¬ 
tally demonstrated in Ref. Hfl 

In the non-resonant (dispersive) regime, the energy shift 
of the super-radiant mode is of order G 2 / A, where G is 
the collective SE-cavity coupling and A = w_i — w c (0) 
is the detuning between the resonator frequency and the 
spin gap. By assuming a spin ensemble with gaussian 
broadening and standard deviation cr, the cavity protec¬ 
tion condition is fulfilled if G 2 /A a. However, reduc¬ 
ing the detuning (or, equivalently, increasing the coupling 
strength) leads to unwanted oscillations of a significant 
fraction (~ G/ A) of the wave-function between logical 
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FIG. 4: Quantum circuit description for the simulation of the hopping part of the Fermi-Hubbard model on a two-dimensional 
lattice. Here we explicitly show the case of i-L\ = —A (c[cs + CgCi^, with ip = 2Ar. R x (9) and R y (6) indicate single-qubit 
rotations about x or y axis of the Bloch sphere, while <F((p) is the single-qubit phase gate. 


states |0) and |1). A trade-off can be found in the limit 
of very large G (200-300 MHz) and A (~ 15 — 20 GHz), 
but this is not within reach of present technology. Nev¬ 
ertheless, we show that these oscillations do not prevent 
to implement the proposed digital simulation scheme, 
thus making it possible to employ experimentally avail¬ 
able conditions. In the following we use small detunings, 
i.e. A = 6(5, with G = 30 MHz, and a SE with gaus- 
sian broadening and FWHM = 1 MHz, as assumed in 
S-'I'tI . We numerically determine the time evolution of the 
qubit wave-function, coupled to a bath of dark modes, 
by exploiting the formalism developed in [37;]. We ob¬ 
tain a wave-function leakage to dark states at long times 
of only ~ 1 %. As expected, this is lower than the up¬ 
per bound (4 ct 2 A 2 /G 4 ~ 2.8%) obtained in [38}. As we 
explicitly show below, the wave-function oscillations as¬ 
sociated with the relatively large value of (5/A can be 
easily adjusted in our scheme, since they are coherent 
and correspond to single-qubit rotations. 

Here we proceed to illustrate the use of our simulation 
scheme in a cavity-protected regime and specific SE pa¬ 
rameters. The best spin systems are the so-called S'-ions 
(such as Fe 3+ or Gd 3+ ) whose orbital angular momen¬ 
tum vanishes because of Hund’s rules. This makes the 
ion practically insensitive to disorder in the environment. 
In addition, the number of nuclear spins should be mini¬ 
mized, since these produce random quasi-static magnetic 
fields causing IB. Linewidths as small as a fraction of 
MHz are indeed observed in diluted magnetic semicon¬ 
ductors, such as Fe 3+ in ZnS [39[, whose nuclei are mostly 
spinless. Even if Fe 3+ and Gcl 3+ are S > 1 spins, these 
behave as effective S = 1 SEs if all magnetic-dipole tran¬ 
sitions from the ground stated are far detuned from res¬ 
onator frequencies, apart from two. 

To keep the experimental demonstration of the proposed 
scheme as easy as possible, we choose resonator frequen¬ 


cies smaller than in Section II. In particular, we assume 
14 GHz for the logical , and 10.2 GHz for the auxiliary 
resonators, respectively, which are comparable with the 
resonator frequencies already employed in Ref. [23|. Spin 
ensembles fitting our scheme with a 14 GHz supercon¬ 
ducting resonator can be easily found. For instance, 
Fe 3+ impurities in the same AI 2 O 3 matrix employed 
in [23| display suitable gaps with an applied magnetic 
field of ~ 70 mT forming an angle of ~ 70° with the 
anisotropy axis (given an easy plane axial anisotropy, 
with D = 5.15 GHz 0 )- As mentioned above, the rel¬ 
atively small detuning A = 6(5 needed for cavity pro¬ 
tection induces an unwanted one-qubit oscillation with 
frequency v = \/G 2 + A 2 /4 (~ 95 MHz), which must 
be taken into account in our simulation scheme. This 
can be done in the implementation of one-qubit gates by 
choosing a starting time of the gate ti = 2n'K/v and by 
making the detuning A G for the short time (< 1 ns) 
needed to implement the phase gate, Eq. IB1I We stress 
that this time is much shorter than that characterizing 
the damping due to IB. As far as two-qubit gates are con¬ 
cerned, the only part which is affected by the unwanted 
oscillations is the photon hopping process between logi¬ 
cal and auxiliary cavities, whose starting time needs to be 
chosen again as ti = 2mr /V (note that the auxiliary res¬ 
onator does not contain a SE). In some situations (such 
as the simulation of T-L xz terms), when a pair of qubits is 
subject to a different or differently ordered sequence of 
operations, we need to increase the detuning for one of 
the two qubits to freeze its evolution, while waiting the 
other to complete an oscillation (”rephasing”). 

To test the performance of this setup we first numeri¬ 
cally determine the fidelity of simulations of the elemen¬ 
tary terms in a one- and two-qubit Hamiltonian, includ¬ 
ing also decoherence effects (by solving the Liouville-von 
Neumann equations of motion). Table Q] shows that these 
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FIG. 5: Simulation of the XY model on a pair of hybrid 
qubits using the cavity-protected regime described in the text 
(' Hxy = A[six S2x +si w S2y])- Lines represent the exact evolu¬ 
tion, whereas points are calculated with Lindblad formalism 
assuming Q = 10 6 and for two different values of Tff . We use 
the parameters cu c (0)/27r = 14 GHz, cj_i/27t = 14.18 GHz, 
cji/27r = 12 GHz, w c _(0)/2tt = 10.2 GHz, fi 0 i/27r = 9.2 GHz, 
12 12 /27r = 8.3 GHz. G_i = 30 MHz, G+i = 33 MHz, A = 6G. 


fidelities remain high, thus demonstrating the effective¬ 
ness of our scheme even in a cavity-protected regime. An 
interesting proof-of-principle experiment that could be 
readily performed with existing setups is the simulation 
of the dynamics resulting from an XY interaction be¬ 
tween two spins s = 1/2. This is also the central step 
in the simulation of hopping processes in fermion Hamil¬ 
tonians (dashed boxes in Fig. HJ). Fig. [5] shows that 
the time evolution is well reproduced even with a short 
r 2 tr = l pa. 


V. DISCUSSION AND PERSPECTIVES 

We have proposed a digital quantum simulator based 
on hybrid spin-photon qubits, encoded in an array of 
superconducting resonators strongly coupled to spin en¬ 
sembles. Within this quantum computing architecture, 
quantum gates are implemented by a single operational 
tool, namely by tuning the resonators frequencies. We 
have shown the feasibility of the scheme with state-of 
the-art superconducting arrays technology, which allows 
the high fidelity simulation of a large class of multi¬ 
qubits spin and fermionic models. To test our predic¬ 
tions, we have performed numerical simulations of the 
master equation for the system density matrix, including 
the most important decoherence channels such as photon 
loss and pure-dephasing of the transmon involved in two- 
qubit entangling gates. Coherence times of single spins 
are so long that their effect on quantum simulations can 
be disregarded. However, inhomogeneous broadening of 
the SE might result in an irreversible population leakage 
out of the computational basis. Nevertheless, we have 
shown in the previous Section that this problem can be 


circumvented by implementing the scheme in a cavity- 
protected regime. 

Sources of errors. We analyze here the sources of error 
that affect the quantum simulation, and point out possi¬ 
ble solutions. Three main simulation errors can be found: 
digital errors (arising from the Trotter-Suzuki approxi¬ 
mation), gating errors (due to imperfect implementation 
of the desired unitaries), and decoherence errors (due to 
the interaction of the quantum simulator with the envi¬ 
ronment). While digital errors can obviously be reduced 
by increasing the number of Trotter steps, gating errors 
are accumulated by repeating a large number of quantum 
operations. In the same way the effects of the interac¬ 
tion of the system with the environment become much 
more pronounced if the number of gates, and thus the 
overall simulation time, increases. First, we notice that 
the present setup limits the role of the transmon , which 
is not involved in the definition of the qubits. All trans- 
mons are kept in their ground states apart for the specific 
transmons involved in two-qubit gates, which are excited 
only for a short time. Thus, typical state-of-the-art tech¬ 
nology, which ensures transmon dephasing times of the 
order of tens of microseconds, is sufficient to obtain high 
fidelity quantum simulations of relatively large systems. 
Indeed, the three-spins TIM model reported here can im¬ 
mediately be extended to simulate spin Hamiltonians in¬ 
volving many spins, by addressing the different cavities 
in parallel. Photon loss represents the main source of 
decoherence in our hybrid dual-rail encoding. Finally, 
gating errors are mainly due to the small difference be¬ 
tween the photon frequency and transmonic gaps in the 
auxiliary cavities, which induces a residual interaction 
that is never completely switched off. Here we use the 
tunability of the resonator frequency as the only tool to 
process quantum information, but the flux control of the 
Josephson energy of the transmons can also be exploited 
to increase the detuning, thus leading to even larger fi¬ 
delities. 

Two-dimensional arrays. While any model can be im¬ 
plemented onto a one-dimensional register (e.g., the one 
schematically illustrated in Fig. [TJ at the cost of requir¬ 
ing long-range two-qubit gates, it is clear that a regis¬ 
ter topology directly mimicking the target Hamiltonian 
would greatly reduce the simulation effort. In particu¬ 
lar, there are several important Hamiltonians defined on 
two-dimensional lattices whose simulation would greatly 
benefit from a two-dimensional register. Here, we point 
out that our scheme is straightforwardly usable on such 
a register, but its experimental realization necessarily re¬ 
quires the implementation of two sub-lattices of cavities, 
alternatively coupled to spin and transmon qubits, re¬ 
spectively. Fortunately, resonator arrays with complex 
network topologies are realistically possible, already, as 
each cavity can easily couple to multiple other resonators. 
Fig. [6] displays the schematic drawing of a potential two- 
dimensional layout showing how such sub-lattices could 
feasibly realize a two-dimensional simulator. ^From a 
technological point of view, we notice that similar lat- 
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FIG. 6: (Color online). Schematic representation of a two- 
dimensional implementation of the digital quantum simula¬ 
tor. Dark lines show superconducting coplanar resonators 
routed such that each resonator is coupled to four adjacent 
resonators. Yellow boxes indicate logical resonators contain¬ 
ing ensembles of S=1 spins near the magnetic field antin¬ 
odes, while green boxes indicate auxiliary resonators contain¬ 
ing transmons near voltage antinodes. Flux biasing of the 
resonator SQUIDs could be accomplished using microwave 
lines placed on another layer. 


tices with transmon qubits have been fabricated with 
more than 200 coupled cavities. While local tuning in 
such a lattice would require local flux bias on a sepa¬ 
rate layer, this need for local control lines applies to any 
adjustable quantum simulator. On the other hand, we 
notice that a recent technology has shown promising re¬ 
sults to bring flux lines to the interior part of a lattice 
made of a small number of nodes, e.g. by using Alu¬ 
minum airbridge crossovers to route microwave signals 
into a target resonator [4lj . 

Summary. In conclusion, the proposed setup exploits 
the best characteristics of distinct physical systems: the 
long coherence times of the spins, which can encode quan¬ 
tum information and protect it from decoherence, and 
the mobility of photons entering this hybrid encoding 
of qubits. In the end, this allows to realize long-range 
two-body interactions between distant qubits without the 
need for much more demanding SWAP gates. Moreover, 
on-site tunability and scalability make this architecture 
extremely appealing and competitive with respect to al¬ 
ternative proposals, either based on superconducting ar¬ 
rays or on different technologies. 
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Appendix A: Hybrid dual-rail encoding 

We consider a coplanar waveguide resonator containing 
a single photon in a mode of frequency lo c , and an ensem¬ 
ble of A f non-interacting and equally oriented s = 1 spins. 
In the low-excitation regime, the SE can be modeled by 
two independent harmonic oscillators, related to two dif¬ 
ferent magnetic-dipole transitions from the m = 0 ground 
state of the single spin, to the m = — 1 and m = 1 states, 
with excitation frequencies w_i and uq. This can be 
achieved by properly choosing a system with easy-plane 
magnetic anisotropy, which provides a zero-field splitting 
between the m = 0 ground state and the excited m = ±1 
doublet, and in the presence of a small static magnetic 
field. We suppose to initialize the system by preparing 
each spin in its ground state: |0o) = |0i...0jy)- 
If the resonator frequency is tuned to match the spin 
gap uq, the SE can absorb the photon and collectively 
evolve into the state |0i) = |0i...l g ...0jv). 

Transitions between |0 q) and are described (in the 
limit of low number of excitations) by the bosonic op¬ 
erators Si and SI, where Si = |0)(1|<? and 

[Si,si] = i nam. Conversely, if the resonator fre¬ 
quency is tuned to w_i, the SE can evolve into the 
state |0_i) = |0i... — 1 9 ...0 aa), the transition 

being described by the operators S_i and S_ 1; where 

Within the single-excitation subspace of the system 
formed by the cavity mode and the SE, we introduce 
the hybrid dual-rail encoding of the qubit /i: 

|0) M = &-i, At |0) = l^-i>^ = O), 

|l)^at|0) = |^,n M = l), (Al) 

where at is the photon creation operator and |0) = 
|0o, = 0) is the vacuum state. 

Appendix B: Single- and two-qubit gates 
1. Single-qubit rotations 

Resonant processes involving the absorption (emission) 
of the photons entering the hybrid encoding in (Eq. IA1I) 
are exploited to perform one- and two-qubit gates. These 
processes are induced by “shift pulses”, in which the fre¬ 
quency of cavity /i is varied by a quantity for a suitable 
amount of time. In the idle mode, the photon frequencies 
are largely detuned from the spin energy gaps, and Hint 
is ineffective. In addition, the modes and Co 3 c of neigh¬ 
boring cavities are far-detuned and the effect of H p h- P h 
is negligible. Single-qubit gates can thus be performed 
independently on each qubit, which can be individually 
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addressed. 

Off-resonance pulses are employed to obtain a rotation by 
an arbitrary angle about the z axis of the Bloch sphere. 
These induce a phase difference between the |0) and |1) 
states of the hybrid qubits (Eq. IA1D and performs the 
well-known phase gate: 


mT) = ( 0 e-2fr ) ■ (B1) 

where we have assumed step-like pulses of amplitude 
and duration T. 

Conversely, resonant pulses are employed to transfer the 
excitation between SEs and resonators. This produces a 
generic rotation in the x-y plane of the Bloch sphere: 

D (fj\ = ( cos (9/2) -i e~ iS c to sin (9/ 2)\ 
xy sin (0/2) cos (0/2) y 

( B2 ) 

with 0 = G_iT. By properly tuning the initial time 
we can obtain rotations about x (Sgto = 2kn) or y 
(5gt o = (4 k + 1)7t/ 2) axis, while t he p ulse duration 
controls the rotation angle. See Ref. [lf| for a detailed 
derivation. 


2. Controlled-phase gate 

The Controlled-phase (Cp) two-qubit gate is repre¬ 
sented by the matrix: 


/ 

1 

0 

0 

0 ^ 


0 

1 

0 

0 


0 

0 

1 

0 

l 

0 

0 

0 

e~ iv / 


(B3) 


It can be implemented by means of two-step semi¬ 
resonant Rabi oscillations of the transmon state between 
\ipo,j) and \i/) 2 ,j)- We describe here the C p multi-step 
pulse sequence on two qubits initialized in the state 
|l M lj,), as schematically shown in Fig. [U-(b) for /i = 2, 
v = 3 and j = 2: 

1. The first step corresponds to the hopping of the 
photon from logical cavity 3 to the auxiliary res¬ 
onator 2 (interposed between qubits 2 and 3), by 
means of a 7r-pulse that brings the two cavities into 
resonance. 


2. As a second step, the frequency of resonator p = 2 
(w/) is tuned to floi by means of a 7r-pulse, which 
transfers the excitation to the intermediate level 
\tpij= 2 ) of the transmon. 

3. A 7r-pulse is exploited to induce the hopping of a 
second photon from logical cavity 2 to the auxiliary 
resonator. 


4. Then, a semi-resonant process (during which the 
resonator is detuned from the transmon gap by 
a small amount J 12 ) is exploited to induce an 
arbitrary phase on the 11213 ) component of the 
wavefunction [li|. A pulse of duration At = 
1 . 77 ., where <$12 = fli 2 — A/, is the detun- 

V Gl2+‘ 5 12/ 4 

ing between the resonator mode and the \ipip) —> 
lf/ 2 , 2 ) transition of the transmon, adds a phase 
to the system wavefunction. 


<p = 7T — 7T 


S 12 


V | 5?2+ 4G i: 


5. Finally, the repetition of the first three steps brings 
the state back to 11213 ), with an overall phase p. 
By properly setting the delay between the two it 
pulses corresponding to the previous steps (or by 
performing single-qubit phase shifts), the associ¬ 
ated absorption and emission processes yield a zero 
additional phase. 


Conversely, the other basis states do not acquire any 
phase, as required for the Cp gate, due to the absence 
of at least one of the two photons (see Ref. El)- For 
<5i2 = 0 we obtain the usual full Rabi process, which im¬ 
plements a Controlled-Z (CZ). 

The setup is simplified with respect to our previous pro¬ 
posal |15j . as each resonator contains a single photonic 
mode. 

It is also important to note that here we are using an 
ensemble of effective spins S' = 1 as this ensures the pos¬ 
sibility of implementing Controlled-phase gates between 
distant qubits, with no need of performing highly de¬ 
manding and error-prone sequences of two-qubit SWAP 
gates. Long-distance two-qubit interactions are a key- 
resource for the digital simulation of many interesting 
physical Hamiltonians. They appear each time that a 
multi-dimensional target system is mapped onto a lin¬ 
ear chain of qubits or in models with Wbody terms. 
Among these, as discussed in Section III.B, a particu¬ 
lar interest is assumed by problems involving interacting 
fermions in two or higher spatial dimensions, which are 
often intractable for classical computers. For instance, 
solving the two-dimensional Hubbard model is considered 
by many as the ultimate goal of the theory of strongly 
correlated systems. In these cases the Jordan-Wigner 
mapping induces many-spin interactions pfj ] which can 
be handled as outlined in Fig. 2J provided the ability to 
efficiently implementing long-range two-qubit couplings. 
These are obtained by bringing the photon components 
of the two qubits into neighboring logical resonators by 
a series of hoppings. The operations outlined in Fig. |T} 
(b) are then performed to implement a Cp gate between 
neighboring qubits, and the photon components are fi¬ 
nally brought back to the starting position by reverting 
the series of hoppings. The photons can be transferred 
with negligible leakage and without perturbing the inter¬ 
posed qubits by temporarily storing the photon compo¬ 
nent of these qubits into the m = 1 spin oscillator. We 
stress that a large number of these long-range two-qubit 
gates can be implemented in parallel in the actual setup. 
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Appendix C: Density matrix master equation 


The time evolution of the system density matrix p 
is described within a Markovian approximation and a 
Lindblad-type dynamics, with the Liouville-von Neu¬ 
mann equation of motion [42j : 




+ E r A[p]+E^( C1 ) 

i i 


being r ? - and 7 j respectively the damping and pure- 
dephasing rates of the field Xj. The Lindblad term for 
an arbitrary operator, x, is given by 

Cx\p\ = ~ 2 (x^xp + px^x) + xpx^. (C2) 

If the operator Xi destroys an excitation in the sys¬ 
tem, terms like vCaJp] account for energy losses, while 
pure dephasing processes are described by £ .[p]. We 
note fl 6 j that the former ones provide the most impor¬ 
tant contribution for photons 01 (with x t = a M , cij ), 
while the latter are very important for the transmons 
(sbi = |V’fc,j)(V , fc+ijl, & = O^l)- We represent each field 
as a matrix in the Fock-states basis, and truncate it at a 
number of total excitations previously checked for conver¬ 
gence. The total Hamiltonian, Eq. m, and the density 
matrix master equation of the whole system, Eq. m, 
are built by tensor products of these operators. Then, 
the equation of motion for p is numerically integrated, in 
the interaction picture, by using a standard Runge-Kutta 
approximation. 


Appendix D: Interacting spin fermions 

To extend the quantum simulation of two-dimensional 
Hubbard models to the case of fermionic systems with 
spin, we need to encode each fermion operator into a 
pair of qubits, corresponding to spin up and spin down. 
To achieve this, we exploit a generalization of the Jordan- 
Wigner transformation [44| . For this mapping we need to 

introduce two different spin 1/2 operators S and T, with 
SZ»-i = cj‘ t c Mt - \ and = c^c m4 . - |, describing 
respectively odd and even qubits (ordered by rows in the 
two-dimensional lattice). 

C+ _ Jir T,u=i s ij. 6 ^+* 7r £2=1 4 

(Dl) 


It can be shown that these operators satisfy the 
usual angular momentum commutator algebra, and that 




= 0. We assume that the fermion variables 


are ordered by rows in the Hamiltonian. The efficiency of 
the scheme would be increased by using a 2 -dimensional 
setup consisting of N rows and 2 M columns. We can 
write the Hubbard Hamiltonian in terms of the spin vari¬ 
ables introduced above 


NM 

H s Hub =-A E (-1 )« + 1 (^- 1 +^)5+_ i 5-_ i + h.c. 

NM 

-A E (_i)E^ + 1 (^)f+f 2 -+h.c. (D2) 


NM 


+uJ2 +| E (4Vi +%) 


NM 


+ - 


NMU 


11=1 


M=1 


4 ’ 


where p and v are nearest neighbors on the two- 
dimensional fermionic lattice, such that v = p + 1 (hori¬ 
zontal neighbors) or v = p + M (vertical neighbors) with 
the present labeling. Odd (even) qubits encode spin up 
(spin down) variables. Since the hopping term does not 
act if c^v) = 1 (i.e. cjj^ = 0 ), we can start directly 
with 7 = p + 1 , and the exponential in expressions like 
5+ exp{i 7 r (£7 + ^)}E can be factorized. We 
note that in the case of horizontal neighbors the phase 
factor cancels out and that in Hfj ub do not appear terms 
Stu-iTzv, as we are n °t considering spin-flip processes. 
To simulate such evolution we can proceed in a way anal¬ 
ogous to the spinless case. Here, however, two different 
series of CZ Mi7 should be carried out, depending if we 
are considering the hopping of spin f or spin 4 * fermions. 
The former involves only odd values of 7 , the second only 
even. Notice that, in a 2-dimensional register, we need 
to transfer photons to implement or T^T^, 

each time we have to couple a pair of fermions belonging 
to the same row (due to the alternating f-| mapping), 
but in that case ]/[ Ucz,, 7 is not required. The term 
[\UcZu.~ j needed to correct the sign problem , is neces¬ 
sary only if v = p + M (no photon transfer in that case 
is needed). 
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